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Abstract A quasi-non-linear fracture mechanics model 
based on beam on elastic foundation theory is applied for 
analysis of dowel joints with a single dowel loaded perpen¬ 
dicular to grain. The properties of the elastic foundation are 
chosen so that the perpendicular-to-grain tensile strength 
and fracture energy properties of the wood are correctly 
represented. It is shown that this particular choice of foun¬ 
dation stiffness makes a conventional maximum stress fail¬ 
ure criterion lead to the same solution as the compliance 
method of fracture mechanics. Results of linear elastic frac¬ 
ture mechanics are obtained as a special case by assuming 
an infinitely large value of the foundation stiffness. Results 
of tests on so-called plate joints are compared with theoreti¬ 
cal predictions, showing good agreement for variations in 
initial crack length as well as edge distance. 
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Introduction 

Dowel-type fastener joints subjected to loading perpendicu¬ 
lar to grain may fail in either a ductile manner, character¬ 
ized by bending of the fastener and/or embedment of the 
fastener into the wood, or in a brittle manner characterized 
by splitting of the wood. The ductile failure modes are well 
understood and can be predicted with good accuracy by the 
European Yield Model, originally proposed by Johansen, 1 
which now forms the basis of design of dowel-type fastener 
joints in major design codes. The brittle splitting failure 
mode is, however, not yet well understood. Although a few 
simple models suitable for implementation in design codes 
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have been proposed recently, 2-7 none of those models have 
yet gained wide acceptance. 

Plate joint specimens as shown in Fig. 1 have previously 
been used to derive relevant information on dowel-type 
fastener joints loaded perpendicular to grain. 5-6,8 The model 
presented in this article considers a double symmetrical 
plate joint specimen with a single dowel, but may also be 
applied to joints in beams for which the edge distance, h c , is 
small as compared with the total beam depth. The model is 
based on beam on elastic foundation (BEF) theory for a 
Timoshenko beam, and the foundation stiffness is asso¬ 
ciated with the fracture performance of the wood, i.e., 
perpendicular-to-grain tensile strength and mode I fracture 
energy. 

A model based on the same ideas as the present model 
was previously reported 9 considering joints with multiple 
dowels in a row parallel to grain. The previous model in¬ 
cludes effects of edge distances, end distances, and dowel 
spacing. However, it does not address the effect of crack 
length, including dowel hole size, initial drying cracks, 
cracks occurring due to local compression damage in the 
wood beneath the dowel, etc. The effect of crack length is 
the main point of focus in the present article, which is lim¬ 
ited to considering joints with a single dowel and with infi¬ 
nite end distances. Furthermore, the previous model 9 was 
derived using conventional stress analysis and by introduc¬ 
ing fracture mechanics through the constitutive relation of 
the foundation. In the present article, full derivations are 
given based on the compliance method of fracture mechan¬ 
ics and on the conventional stress approach, and it is shown 
that a certain choice of the foundation stiffness leads to the 
same predicted failure load with either approach. 


Theory 

The symmetrical half of the upper joint in Fig. 1 is modeled 
as a beam on elastic foundation as schematically shown in 
Fig. 2. A state, in which a crack of length 2 a has developed, 
is considered. Let w 0 = w(0) and 0 O = 0(0) denote the 
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Fig. 1. Geometry of plate joint specimen 
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Energy release rate failure criterion 

The deflection at the loading point, <5 P , may be written 
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and d Pw and are the contributions from the beam on 
elastic foundation 
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The compliance, C, is given by C = d P IP, and the increase in 
compliance due to an increase in crack length is found to be 
given by 
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Fig. 2. Symmetrical half of plate joint modeled as beam on 
elastic foundation 


displacement and rotation of the beam axis at x = 0, respec¬ 
tively, with positive directions as indicated in Fig. 2. 

From the theory of a Timoshenko beam on elastic 
foundation 1011 
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For a linear elastic body loaded by a single load, P, the crack 
propagation energy release rate, Gj. is given by 12 


P^dC 

* 2b dfl 
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A crack starts propagating when the energy release rate 
assumes a critical value, (g„ i.e., the failure criterion is 


5= & 


( 10 ) 


and where P 0 and M 0 are shear force and moment in the 
beam at x = 0, respectively, b is beam width, E is the 
modulus of elasticity (MOE) in the grain direction of 
the beam, G is the shear modulus of the beam, K is the 
stiffness of the elastic foundation (units: N/m 3 ), I is moment 
of inertia, and A is cross-sectional area of the beam (rectan¬ 
gular cross section is assumed). 

For the structure shown in Fig. 2, P 0 and M 0 are found to 
be 


Assuming static or quasi-static conditions and no energy 
dissipation outside the fracture region, the critical energy 
release rate is equal to the material property fracture en¬ 
ergy, g f , i.e., 


(3) 


From Eqs. 9-11 it follows 
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Especially for a —» 0, Eqs. 8 and 12 lead to 
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Maximum stress failure criterion 


The maximum tensile stress in the foundation occurs at 
x — 0 and is given by 

o Q = Kw 0 (14) 

Failure is assumed to occur when the maximum tensile 
stress equals the perpendicular-to-grain tensile strength, 
i.e., 

0o=/ t (15) 

From Eqs. 1, 3, 14, and 15 it follows that 


2 bft p + (1 — K )X 2 a 
Especially for a —» 0, Eq. 16 leads to 

P - = bft 2 /- l 1 


(16) 


(17) 


Foundation properties 


The deformation of the foundation is, in general, assumed 
to be composed of a contribution from a special fracture 
layer (with no physical thickness) modeling the fracture 
damage and a contribution from the perpendicular-to-grain 
elastic strains in the wood as indicated in Fig. 3. 

The damage and fracture performance of wood is in 
general nonlinear, but is in the present analysis represented 
by a linear response that is equivalent in terms of peak 
stress, f, and fracture energy dissipation, (g t . Because the 
special fracture layer has no physical perpendicular-to-grain 
dimension, strain is not defined and the perpendicular-to- 
grain stress, o, is used as constitutive relation as a function 
of the perpendicular-to-grain deformation, d f . 

From Fig. 4, it follows that the stiffness of the fracture 
layer, K f , is given by 


K t = 

2<3i 


(18) 


axis of beam on elastic foundation 



Fig. 3. Contributions to foundation stiffness 



The elastic perpendicular-to-grain strain, e s , in the wood 
is given by Flook’s law, and the deformation, d s , and stiff¬ 
ness, K s , of the part of the specimen between the beam axis 
and horizontal line of symmetry is given by 
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The slice of the foundation shown in Fig. 3 is assumed to be 
in a state of uniaxial stress, and the total stiffness of the 
foundation is thus given by 


K = 


K t K s 
K t + K s 


( 20 ) 


In previous applications of beam on elastic foundation 
theory to mode I fracture problems, 11,1314 the special frac- 


Fig. 4. Constitutive relation of fracture layer 


ture layer as considered in the present article was omitted, 
and the deformations of the foundation have been attrib¬ 
uted solely to the elastic perpendicular-to-grain strains. If 
disregarding the elastic perpendicular-to-grain strains, i.e., 
K s —» °°, then 



( 21 ) 


Using the stress approach as given by Eq. 16 or 17 
together with Eq. 21 is a complete analogy to the fracture 
mechanics application of the Volkersen model to analysis of 
mode II failure in lap joints. 15 The analysis has been termed 
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quasi-non-linear fracture mechanics because the material 
responses are assumed to be linear as in linear elastic frac¬ 
ture mechanics (LEFM). However, at the same time, the 
tensile strength is assigned a finite value, not an infinite 
value as in LEFM, and, while it is a characteristic of LEFM 
that the failure load is proportional to the square root of the 
fracture energy, this is not a characteristic of the quasi-non- 
linear fracture mechanics model. Equation 21 ensures that 
the fracture mechanics and stress approaches lead to the 
same failure load. In general, when using BEF theory it is 
significantly easier to derive simple solutions by means of 
the stress approach than by the fracture mechanics ap¬ 
proach. Equation 21 makes Eqs. 13 and 17 lead to the same 
solution as previously derived for glued-in rods subjected to 
pure shear 9 ' 16 and for the splitting strength of dowel joints 
loaded perpendicular to grain. 16 This solution may be writ¬ 
ten (note that the total failure load of the plate joint speci¬ 
men according to Fig. 1 is 2 P c0 ). 


Figure 6 shows y = /’ d // J I l;l kl as given by Eq. 22 as a 
function of the perpendicular-to-grain tensile strength. The 
same b, E, G, and Cj, as used in Fig. 5 apply to Fig. 6. 
The sensitivity to variations in/ t is similar for finite values of 
the initial crack length (Eqs. 16 and 21). 


Experimental 

Three test series were conducted on plate joint specimens. 
Series 1 was conducted in order to evaluate the influence of 
the initial crack length. Series 2 was conducted to evaluate 
the influence of the edge distance. Series 3 was conducted 
to evaluate the influence of elastic perpendicular-to-grain 
strains. 

Methods and materials 
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For f or E —» °°, the failure load according to Eq. 22 be¬ 
comes 2 P c0 = 2h(5GGjfliJ3y 2 . This is the same solution as 
previously obtained by using a fracture energy balance 
equation and a chosen displacement field, which considers 
only shear deformations, 5 ’ 6 if in the latter solution the shear 
area A s = 5bhJ6 is used instead of the full cross sectional 
area A = bh e , and the same solution as predicted for hJh—> 0 
by another previous LEFM solution. 2,7,8 

Figure 5 shows an example of the influence of crack 
length, a, on the failure load according to Eqs. 8 and 12 
or Eq. 16. Equation 21 has been assumed, and the mate¬ 
rial and geometrical properties used are: h s = 40mm, 
b = 25 mm, E = 5670MPa, G = 315 MPa, f t = 3.5 MPa, 
( g f = 0.20 N/mm. 
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Fig. 5. Influence of crack length on the failure load 


All specimens were made of glulam of Japanese cedar 
(Cryptomeria japonica ). The MOE in the grain direction 
was determined by measuring the longitudinal vibration 
frequency of the glulam beams, from which the specimens 
were cut. All tests were displacement controlled and time to 
failure was 2-3 min. All tests were conducted using a 14-mm 
dowel in a 15-mm hole. Specimens were not especially se¬ 
lected to avoid knots or other defects. 

Series 1 was conducted to evaluate the influence of the 
initial crack length, a. All specimens had dimensions of b = 
25 mm, h = 200 mm, h c = 40 mm, and L — 500 mm. Ratios 
alh s = 0, 0.5, 1, and 2 were tested. The moisture content 
(MC) was 10.5%, density was 370kg/m 3 at 10.5% MC, and 
the MOE was 5670 MPa. Twenty specimens were tested for 
a//i e = 0, and ten specimens were tested for other a!h c ratios. 
The cracks were cut with a small band saw resulting in a 
crack width of 1 mm. The reported crack length is the length 
of the cut crack, i.e., from crack tip to the dowel hole 
periphery. Tensile strength perpendicular to grain was 
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Fig. 6. Normalized failure load as a function of perpendicular-to-grain 
tensile strength for zero crack length 
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Fig. 7. Test specimen for determination of perpendicular-to-grain ten¬ 
sile strength 



Units: mm 


Fig. 8. Double cantilever beam specimen for determination of fracture 
energy 


determined using specimens as shown in Fig. 7. Twenty- 
four specimens were tested. Fracture energy was deter¬ 
mined using DCB specimens as shown in Fig. 8, and using 
Eq. 23 with b = 25 mm, h = 40 mm, a = 200 mm, and G = 
E! 18. 
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Series 2 was conducted to evaluate the influence of edge 
distance, h e . Specimens were cut from the same beams as 
the specimens used in series 1. In series 2, no initial cracks 
were cut, i.e., a = 0. All other conditions were the same as 
for series 1. Twenty specimens were tested for h c = 40 mm 
(from series 1), and ten specimens were tested for h e = 
20 mm and for h e = 60 mm. 

Series 3 was conducted to evaluate the influence of the 
elastic strains in the wood. All specimens had dimensions of 
b = 40 mm, h c = 56 mm, L = 500 mm, and no initial cracks 
were cut, i.e., a = 0. Two values of specimen depth were 
tested: h = 800 mm (seven specimens) and h = 200 mm 
(ten specimens). The MC was 12%, and the MOE was 
11363 MPa. 


Results and discussion 

Series 1 



a/h e 

Fig. 9. Calculated and tested failure loads as a function of the initial 
crack length 



h e [mm] 

Fig. 10. Calculated and tested failure loads as a function of the edge 
distance 


ized failure load, PJP B0 , as a function of the normalized 
initial crack length, alh e . P c0 for the tested values is taken as 
the mean value of failure load for a = 0mm ( P d) = 3.15 kN). 
Note that a here denotes the length of the (half) crack 
cut by saw, i.e., the dowel hole is not counted as a part of 
the initial crack. If the dowel hole is counted as part of the 
initial crack, the tested alh e values thus become 0.69, 
1.19, and 2.19 instead of the values 0.5,1.0, and 2.0, respec¬ 
tively, as given in Fig. 9; this leads to even better agreement 
between theory and experiment than indicated in Fig. 9. 
The theoretical failure loads have been calculated using 
Eqs. 16 and 17 with the foundation stiffness as given by 
Eq. 21. 


The tension tests gave a value of the perpendicular-to-grain 
tensile strength,/,, of 3.5 ± 1.1 MPa (mean value ± standard 
deviation), and the fracture energy tests gave a fracture 
energy, <g, of 0.20 ± 0.06 N/mm. Figure 9 shows the normal- 


Series 2 

Figure 10 shows the failure load, P c , as a function of the 
edge distance, h c . The theoretical values have been calcu- 
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lated using Eqs. 16 and 21 with/, = 3.5MPa and (g t = 0.20N/ 
mm (as determined by means of tensile and DCB specimen 
test). The maximum embedment at failure observed in 
series 1 and 2 was about 2 mm. 

The perpendicular-to-grain tensile strength of wood de¬ 
pends heavily on the volume subjected to maximum stress, 
and an appropriate volume (or cross-sectional area) of the 
tensile strength specimen is not obvious. It has previously 
been suggested to use the plate joint specimen as a standard 
test specimen for deriving a so-called fracture parameter, 8 
C 1 = (5G@ f /3) 1/2 . Similarly, the foundation stiffness or the 
perpendicular-to-grain tensile strength may, in principle, be 
derived from standard plate joint tests if the fracture energy 
is known. 

From Eqs. 17 and 21 it follows 


K = 


25 G 2 


3 Eh 
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and 
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3 E l 2 
5 b 2 GQA 
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/, = (25) 

As seen from Fig. 6, the predicted failure load asympto¬ 
tically approaches an upper limit for/, —» °° (and thus K —» 
°°). If the predicted failure load is based on the mean value 
of (g t obtained from DCB specimen tests, some individual 
failure loads obtained by plate joint tests may exceed the 
predicted upper limit due to the fact that the fracture en¬ 
ergy of the individual specimen may be considerably higher 
(typically caused by knots) than the mean value. This effect 
leads to £ > 1 and causes numerical problems in Eq. 24. 
Equations 24 and 25 should thus not be used on individual 
plate joint test results, but only on mean values for which 
numerical problems are less likely to occur. 

Using (/, = 0.20N/mm as obtained from the DCB speci¬ 
men tests, Eqs. 24 and 25 lead to / = 3.5 MPa for h e = 
20mm,/, = 3.9MPa for h c = 40mm, and/, = 1.4MPa for h e 
= 60 mm. Figure 6 shows that the quasi-non-linear fracture 
mechanics solution for the plate joint is relatively insensi¬ 
tive to variations in / in the range relevant for glulam. For 
example, the increase in failure load for h e — 20 mm when 
increasing f from 1.4MPa to 3.5 MPa is only about 13%. 
Relatively small variations in failure loads may thus 
cause large variations in the perpendicular-to-grain tensile 
strengths determined by means of Eqs. 24 and 25, 
which explains the significantly lower f value obtained for 
h s = 60 mm. 


Series 3 

P c 200 and P cS00 denote the failure loads of specimens with h = 
200 mm and h = 800 mm, respectively (all other properties 
being the same). The tests gave the following mean values 
± standard deviation: P cloo = 6.16 ± 0.37 kN and P c800 = 6.71 
± 0.30 kN. F-test and Student’s f-test indicate that the two 
specimens can be assumed to have same variation, but dif¬ 
ferent mean values of the failure load. No excessive embed¬ 
ment of the dowel was observed. 


Table 1. Calculated failure load ratios PaJPcsm 


Assumption 

Equation 13 

Equation 17 

k = flag. 

1.00 

1.00 

K = 2 E y l(h - IQ 

1.19 

0.52 


Table 1 shows the calculated failure load ratios, Pciad 
F cg0(1 , given using Eq. 13 (fracture mechanics approach for 
a = 0) and Eq. 17 (stress approach for a = 0) assuming 
either K = f 2 IGj f (i.e., disregarding the elastic strains) or 
K = 2E y l(h ~ h e ) (i.e., disregarding the fracture layer). Note 
that the calculated failure load ratios are all independent of 
/ and <g r . The value of F c200 /F c800 found by testing is 0.92. 
Although the test results of series 3 are too sparse to make 
any firm conclusion, a value of the foundation stiffness inde¬ 
pendent of h as given by Eq. 21 (K = / 2 /2(j £ ) seems to be the 
most appropriate choice. 


Conclusions 

A quasi-non-linear fracture mechanics model for analysis of 
dowel-type fastener joints with a single fastener subjected 
to perpendicular-to-grain loading was presented. The 
model is based on Timoshenko beam on elastic foundation 
theory, which includes shear deformations, and it takes into 
account the crack length. The analysis includes a fracture 
mechanics approach based on the compliance method and a 
conventional stress approach. The properties of the linear 
foundation (fracture layer) are chosen so that the 
perpendicular-to-grain tensile strength and fracture energy 
properties of the wood are correctly represented. This par¬ 
ticular choice of foundation stiffness was shown to unify the 
stress and compliance approaches. Some previous LEFM 
models 257 were shown to appear as special cases of the 
quasi-non-linear fracture mechanics model. The quasi-non- 
linear fracture mechanics model was found to be able to 
predict the influence of edge distance and initial crack 
length in good agreement with test results. 

Quasi-non-linear fracture mechanics models based on 
BEF theory require the perpendicular-to-grain tensile 
strength, /„ as an input, but because f is highly volume 
dependent, direct determination by tension tests requires a 
method for estimation of a proper volume (or cross section) 
of the test specimen. The plate joint specimen has previ¬ 
ously been proposed as a standard test specimen for deriv¬ 
ing relevant fracture properties for use in LEFM models. 5 6 8 
However, the plate joint specimen seems to be too insen¬ 
sitive to variations in the perpendicular-to-grain tensile 
strength. Standard test specimens more sensitive to varia¬ 
tions in/, than the plate joint specimen should therefore be 
sought for determination of the perpendicular-to-grain ten¬ 
sile strength. 
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